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Abstract 

We study the spectrum of density fluctuations in the ultra-cold gas of neutral atoms, confined 
in a magneto-optical trap. We determine the corresponding amplitude and spectra of laser light 
scattered by this medium. We derive an expression for the dynamical structure function, by using 
a test particle method. We propose to use the collective laser scattering as a diagnostic method for 
the microscopic properties of the ultra-cold matter. This will also allow us to check on the atomic 
correlations which are mediated by the collective mean field inside the gas. 
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I. INTRODUCTION 



In recent years, attention has been given to the study of collective oscilations in alkali 
metal gases, cooled and confined in magneto-optical traps (MOTs) [IHS]. Such a cooperative 
behavior is driven by light and results from the random scattering between the atoms, being 
in the root of a long-range interaction potential. Surprisingly, this collective interaction is 
very similar to that observed in Coulomb systems, since it scales as which allows one 
to establish analogies between cold atoms and plasmas. 

The evidence of such Coulomb-like forces due the random scattering of photons was first 
discovered by Walker et al. |1] and Dalibard [S] , where the divergence of the effective net force 
between the atoms was shown to be proportional to the MOT density. As a consequence, 
it is possible to define an equivalent electric charge for the neutral atoms. The advantage 
of such a system is that the effective charge depends on the basic parameters of the system, 
what makes MOTs good candidates to observe plasma physics phenomena with a tunable 
"electric" charge. A strong manifestation of the plasma character of MOTs was also revealed 
by the studies of Coulomb explosions of the atomic clouds, performed by Pruvost et al [B]. 
Recently, we have shown theoretically that waves and oscillations of the plasma type can 
also exist in the neutral gas [3 18] and have explored the plasma character of magneto-optical 
traps to the study of parametric instabilities |9]. 

Here we propose to study the spectrum of density fluctuations in the ultra-cold gas, 
by using collective laser scattering techniques. In doing this, we explore even further the 
analogy between the ultra-cold neutral gas conflned in a magneto-optical trap and a plasma, 
by adapting well known plasma physics theoretical techniques to this new problem. First, 
we determine the scattered flelds, and the expected average scattered power. Second, we 
determine the spectrum of density fluctuations, in the quasi-classical approximation. We use 
here the test particle method, well known in plasma physics (TU]. Finally, we characterize 
the main properties of the expected scattered signal. 

II. SCATTERED LASER POWER 

We start by stating the basic equations of our problem. The propagation of a laser pulse 
along the ultra-cold gas is governed by the wave equation 
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V^-^a5jE=-,„^ (1) 

where E is the laser electric field and P is the polarization field in the gas, as determined 

by 

/■oo 

P(r,i) = eo/ U{r,T)-E{r,t-T)dT, (2) 
where 11 (r, t) is the atomic susceptibility tensor. Using Fourier transformations of the fields 



E(r,c^)exp(-ia;t) ^, (3) 



E 

one obtains 



P(r,a;) = eon(r,a;)-E(r,a;). (4) 

Because we can neglect the anisotropies in the majority of the MOT, the atomic susceptibihty 
tensor 11 can be safely replaced by the spectral susceptibility function % which, for the two- 
level atom, is given by 

Here, n(r) is the number of atoms per unit volume, Xa(i^) is the single atom susceptibility, 
S = ou—Ua is the radiation frequency detuning with respect to the atomic transition frequency 
uja, r is the corresponding line width and fl represents the Rabi frequency. The local density 
can fluctuate in time, at frequencies much lower than those of the radiation field frequency 
u, thus leading to the replacement of n(r) by n{r,t) in the above expression. The total 
radiation field E will be determined by the sum of two parts, the incident laser field Eq, and 
the scattered field E^. For a real incident laser field, we can use 

E(r, t) — ^Eq exp(iko • r — icuot) + c.c, (6) 
where the wave- vector ko follows the dispersion relation 

^0 = ^ [1 + x(^o)] , x(t^o) = noXa(wo), (7) 
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with no being the atomic mean density. Such an approximation is vahd for MOTs with a 
moderate number of atoms (typically, ~ 10®.) On the other hand, the scattered field can 
be generically written as 

/doj 
Es{r,u)exp{-tut)^ (8) 

and the field components can be determined by solving the corresponding wave equation, 
which can be written as 

(V^ + E,(r, tu) = |^Eoe^'^°-'-n(r, to') + c.c, (9) 

where 



n{r,u') = J n{v,t)e"^Ut (10) 

and oj' = oj — ujq is the frequency of density fluctuation. In its turn, the scattered wave 
number k will be determined by the same dispersion relation as in Eq. ([T]), with and ojq 
respectively replaced by k and oj. The solution to Eq. (|9| can therefore be written as 



Es{r,u) = -^^Xa{coo){e^-eo)Eoe^''-^ I dr J ^^nik' ,uy^'^''+^'-^'>-^ + c.c, (11) 

where we have used the unit polarization vectors for the incident and the scattered fields, 
eo = Eq/Eq and e^^ = Es{uj)/ Es{u)). The integration in volume leads to the appearance of 
a delta function 5(ko + k' — k), and subsequent integration over k' finally leads to 



Es{r,ijj) = —iAsexp{ik ■ r) + c.c, (12) 
where the amplitude of the scattered field is determined by 

^2 

Asiuj) = ^^Xa(wo)(e^ • 60)^0^^(0;', k') (13) 

with OJ = uq + u' and k = ko + k'. The relative scattered intensity along the direction of k 
is therefore determined by the dynamical structure factor 

Sik,co) = ^^^1^ = ^IXa(c.o)P(e. ■ eo)Xk',a;')P. (14) 
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This quantity can easily be related with the density correlation function, by noting that 



\n(k',uj')\^ = j dvidti j dv2dt2n{viMy~'^'''''^'"'''' n{v2M)(^^''''^''"'''^- (15) 
Performing an average in both space and time, the latter equation reads 

< |n(k',w')P >= lim — I dridti [ dv2dt2 < n(ri, ti)n(r2, ts) > e^^' e~'^' ^'''-'^^ 
y,T-s>oo VT J J 

(16) 



Equation ( 14 ) establishes the relation between the scattered signal and the fluctuations inside 
the cloud, which contain the basic features of the laser cooled gas. We should therefore 
proceed by understanding how to compute such fluctuations by taking into account the 
microscopic distribution of the atomic states. 



III. TEST PARTICLE METHOD 



In the quasi-classical approximation, a trapped alkali gas in a MOT can be described by 
a kinetic equation of the form [7] 

where W = W{r, v, t) is the atomic distribution function. Here the collective force F can 
be determined by the Poisson equation 

V-F = Qn = Q J Wir, v, t)dv, (18) 

where Q is the atomic effective charge P, [7j 

Q = ^aR{an-aL). (19) 

Here, Iq = 2Jsatf2^/r^ is the total intensity resulting from the six cooling laser beams, I sat 
is the saturation intensity of the atom, c is the speed of light, represents the optical 
cross section for absorbing photons of the laser beams and is the optical cross section for 
absorbing photons that are reradiated from the atoms in the trap. [H]. Our procedure will 
consist in the test particle method [I^ , by considering a single atom inside the trap moving 
with speed v' as the perturbation source. The presence of such a test particle will lead to 
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a fluctuation in the gas distribution function W, as determined by the hnearized kinetic 
equation 



a \ ~ ATT f) 

- + v.v)*-^--.^«^-„. (20) 

where Wq = W — W is the unperturbed (or the equihbrium) distribution function, and the 
collective force perturbation 6F induced by the test particle in the medium is determined 

by 



V-6F = Q j W + 5{r - r')6{w - V) 



(21) 



The test particle position is assumed to be given by r' = r'{t) = fq + v't. Let us multiply 



equation (21) by exp(— ik ■ r + iut) and perfomr integrations over both space and time. 
Deflning the Fourier components 5F^k and P^^k for the perturbed force and distribution, 
respectively, we can easily obtain 



ik ■ (5F^k = 2nQSiLo -k-V) + Q J W^^k(vrfv 
Similarly, the Fourier transformation of equation ( [20| ) leads to 

I 5F^k-9Vro/5v 



M (w-k-v) 



(22) 



(23) 



Replacing this in equation (|22j), and noting that the collective force is purely longitudinal, 
in such a way that (JF^jk = k6Fi^]f^/k, we can write 



5F, 



-1^7, -, k, 



P e(w,k) 

where we have introduced the dielectric function of the atomic gas 



Q/"k - dWo/d-j , 
e{uj, k) = 1 + x(w, k) = 1 + / r^^rfv. 



(24) 



Mfc2 7 (w - k ■ v) 

Replacing equation (|24l) in the expression of the perturbed distribution, we then get 



(25) 



2ttQ 5{uj-k- v') k • dWo/dw 



(26) 



MP e(u;,k) (w - k • v) ' 
After inverse Fourier transforming, we can compute the perturbed force due to the test 
particle 
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dn^.t)- ^ ^2 J e(k-V,k) (27r)3 ^^^^ 

and the associated perturbation to distribution function 

~ Q /• exp[.k-(r-ro-Vt)]k-giyo/gv ^k 
^^"'"^'^^ - ~MF y e(k.v',k) k.(v'-v) (27r)3- ^^^^ 

Finally, the density perturbation created at a position r and instant t by a test particle 

located at position fq at t = moving with velocity v', can be now determined by 

n(r,t) = (5(r-r') + yiy(r,v,t)rfv. (29) 

At this point, it should be notice the power of this method, as all the atoms in the gas can be 
individually considered as test particles. This means that the total averaged perturbations 
can be calculated by integrating the relevant quantities defined above over the equilibrium 
distribution function VFo(ro; v'). In particular, the averaged mean force perturbation can be 
determined by 

<5F{v,t)>= J dvo j dV5F{v,t)Woi^^). (30) 
This quantity is a linear superposition of purely oscillating quantities, and it can easily be 



found using equation (27) that it is identically zero 



< (5F(r,t) >= 0, (31) 

while the averaged quadratic mean force < |5F(r, t)p >, on the other hand, generally is not. 
It will be rather determined by 

2^ f J f J r f dki ki exp[iki • (r — ro — v't)] 



< i*F(M)i^ >^ ^vj " 

dk2 k2 exp[zk2 ■ (r - fq - v't)] 



{2qn)^ki e(k2-v',k2) 

We next express the fluctuations inside the trap in terms of the statistical average over all 
test particles in the system. 
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IV. DENSITY CORRELATIONS 



Let us introduce Wo(v) = rio/o(v), where the new distribution function /o(v) is the 
normahzed distribution function, such that 



j Wo(v)(iv = ^0 y /o(v)(iv = no. 



(33) 



Using equations (28) and (29), it is possible to express the atom density fluctuations due to 



a given test particle initially located at r(t = 0) = fq as 

n(r, t) = 6{r - Fq - Vt) + j W{r, v, t)d\, 
which can be more represented in a more suggestive fashion 



(34) 



n(r, t) 



dk 



exp[ik ■ (r — To — v't)] k • dWo/dv 



(35) 



e(k-v',k) k-(v'-v) 
Here, we have introduced the definition of the equivalent plasma frequency for the neutral 
gas. Up = Quq/M [7j. Using the average as in Eq. (30), we can then state the density 
correlations as 



< n{ri,ti)n{r2,t2) >= uq J dro J dv7o(v')^(ri, ti)n(r2, ^2) 
In explicit form, this can be written in the following way 



(36) 



< n(ri,ti)n(r2,t2) >= '^o / dro / dV fo{v 



dki 



where we have introduced the new auxiliary quantities 



Lfj = ikj ■ (fj - ro - v't) 



and 



dk2 

(M 



:e'^^g{k,,V)g{k2,V) 



(37) 



(38) 



^(kj^V) 



I -tut 



kj ■ dfo/dv 



(39) 



fc|e(k,-v',k,)k,-(v'-v) 
Integration over the initial position for the generic test particle ro leads to the appearance of 
a delta function (5(ki + k2) which, in its turn, allows us to integrate over k2. We then replace 
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ki by k, in order to simplify the notation. By noting that the velocity v only appears in 
the parallel direction with respect to the wavemathbftor k, we can write 



V = M- + V_L 



(40) 



and introducing the parallel distribution function 



we can finally write 



Fo{u) = j /o(v)dv^, 



I j du 1 



(41) 



< n{T^M)n{r2M) >= no / du'Fo{u') / -j^e''^\I{k,u')\\ 



dk 



(2vr) 



(42) 



where we have used the phase 



(p = ik ■ {ri - r2) - iku'{ti - 



(43) 



and new argument 



I{k,u') 



dFo/du 



-du 



(44) 
which 



P J k^e{ku',k){u' -u) 
This is formally very similar to the results obtained for an electron-ion plasma 
is due to the Coulomb like interactions between nearby atoms in the ultra-cold neutral gas 
[7] . The main difference relies on the fact that we only have one particle species here, which 
suggests that the system can be regarded as a one-component plasma. We can now use 



these results to determine the structure factor discussed in Sec. II. Inserting Eqs. (42) to 



(44) in Eq. (16), and successively integrating over ri, r2, ti, ^2, k, and m', we finally obtain 



S{k' , J 



Fo{u:'lk') 



k' 



where we have used 



I{k' , J 



{dFo/du)du 



(45) 



(46) 



P J k'^e{k',u'){u-u'/k') 

In order to illustrate the collective behavior of the system during the scattering, we assume 
that the atoms in the trapped approximately follow a Maxwell-Boltzmann distribution 
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where Vth = ^/^ksT /m represents the thermal velocity of the atoms, fc^ is the Boltzmann 
constant and T is the temperature. Defining the scattering parameter a = k'Xjj, where 
Ad = Vth/ojp is the Debye length, the dynamical structure factor can be explicitly given by 

S{a,u') = — ^ — I{a,u'/up), (48) 
\/7ra 



where a = k'Xo represents the scattering parameter. The latter defines the nature of the 
scattering processes. Moreover, for a < 1, the scattering signal results from a coherent pro- 
cess, as the system is essentially composed of interacting particles. In a plasma language, 
it means that the Debye length A/j is larger or comparable to the perturbation wavelength 
271 /k' and therefore the correlations between the atoms play an important role on the scat- 
tering signal (unscreened, correlated atoms); on the contrary, if a ^ 1, the scattering is 
said to be incoherent and the Debye length is assumed to be much smaller than the typi- 
cal size of the perturbation (screened, uncorrelated atoms). In the figure, we have ploted 



the normalized structure factor in Eq. (48) for three different values of a. We observe a 
strong resonance near u' = Up for the case of coherent scattering {a = 0.85 and a = 1.0), 
while a very broad spectrum is obtained for the case of incoherent scattering (a = 5.0). 
The physical reason for such results are related with the interference mechanism during the 
scattering: the excitation of collective perturbations of the density will result on the con- 
structive interference (and therefore amplification) of the scattered signal around the natural 
frequency characterizing the long-range order of the interaction, i.e. Up. On the contrary, 
if uncorrelated, single-atoms fiuctuations are excited in the systems, the interference will be 
destructive and, therefore, no resonance is expected to be observed in the spectrum S{a, u'). 

Few words are now devoted to discuss the experimental observation of the spectra. The 
scattering mechanism presented in this work is similar to Rayleigh scattering, with the im- 
portant difference that long-range order between the atoms is present, which is naturally 
induced during the laser cooling process. Therefore, this could be experimentally imple- 
mented in magneto-optical traps in order to measure the collective excitations of the system, 
as both coherent and incoherent regimes of scattering are possible, opening the door to a new 
diagnosis. A simple setup for the present technique would consist of a external probe laser 
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FIG. 1: Dynamical structure factor S{a,uj') for different values of the scattering parameter a: blue 
full line (a = 0.85), black dashed line (a = 1.0), and red full line (a = 5.0). Both coherent and 
incoherent scattering processes are represented. 

beam with a very small waist (much smaller than those of the cooling beams) propagating 
inside the cloud. The resulting fluorescent light, modulated by the density fluctuations, can 
be focused into a photodetector placed at a certain angle 9o with respect to the probe beam 
direction. After Fourier transform the photodetector signal, both the scattering parameter 
a and the effective plasma frequency Up can be measured by fitting the fluorescence sig- 



nal with Eq. (48). In order to distinguish the scattering light from the fluorescence signal 
coming from the MOT, the probe beam should operate a large detuning compared to the 
cooling detuning, i.e {\6o\ > \Scooiing\)- 



V. CONCLUSIONS 



In this work we have determined the spectrum of density correlations of ultra-cold atoms 
confined in a magneto-optical trap. We have explored the formal analogies between the 
physics of the neutral atom gas at that of a plasma medium. The existence of such analogies 
is related with the repulsive force between two nearby atoms in the ultra-cold gas, associ- 
ated with the exchange of scattered photons from the laser cooling system. This repulsive 
interaction can be described by Poisson like equations for the collective force, similar to that 
describing electrostatic fields in a real plasma medium. 

Here we have adapted the test particle approach, well known in plasma physics, to the 
case of the neutral cold gas. We have determined the spectrum of density correlations in 
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the gas and we have estabhshed the corresponding spectra of laser hght scattering. We have 
apphed this to a medium in thermal equilibrium and given examples of expected spectra. 
This work proposes the use of laser light scattering as a diagnostic technique for probing the 
internal kinetic properties of a ultra-cold gas, and could be useful for the interpretation of 
future experiments in this field. 
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